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19 Abstract: 

20 The Middle East respiratory syndrome coronavirus (MERS-CoV) non- 

21 structural protein 16 (nsp16) is an S-adenosyl-L-methionine (SAM)- 

22 dependent 2’-0-methyltransferase (MTase) that is thought to methylate 

23 the ribose 2’-OH of the first transcribed nucleotide (Ni) of viral RNA cap 

24 structures. This 2’-0 MTase activity is regulated by nsplO. The 2’-0 

25 methylation prevents virus detection by cell innate immunity mechanisms 

26 and viral translation inhibition by the interferon-stimulated IFIT-1 protein. 

27 To unravel the regulation of nsp10/nsp16 2’-0-MTase activity, we used 

28 purified MERS-CoV nsp16 and nsplO. First, we showed that nsp16 

29 recruited N7-methylated capped RNA and SAM. The SAM binding 

30 promotes then the assembly of the enzymatically active nsp10/nsp16 

31 complex that converted 7m GpppG (cap-0) into 7m GpppG 2 om (cap-1) RNA 

32 by 2’-OFI methylation of Ni in a SAM-dependent manner. The subsequent 

33 release of SAH speeds up nsp10/nsp16 dissociation that stimulates the 

34 reaction turnover. Alanine mutagenesis and RNA binding assays allowed 

35 the identification of the nsp16 residues involved in RNA recognition 

36 forming the RNA binding groove (K46, K170, E203, D133, R38, Y47 and 

37 Y181) and the cap-0 binding site (Y30, Y132 and H174). Finally, we found 

38 that nsp10/nsp16 2’-0-MTase activity is sensitive to known MTase 

39 inhibitors, such as sinefungin and cap analogues. This characterization of 

40 the MERS-CoV 2’-0-MTase is a preliminary step towards the development 

41 of molecules to inhibit cap 2’-0 methylation and to restore the host antiviral 

42 response. 
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43 Importance 

44 MERS-CoV codes for a cap 2’-0-methyltransferase that converts cap-0 into 

45 cap-1 structure in order to prevent virus detection by cell innate immunity 

46 mechanisms. We report the biochemical properties of MERS-CoV 2’0- 

47 methyltransferase, which is stimulated by nspIO acting as an allosteric 

48 activator of the nsp16 2’-0-methyltransferase possibly through enhanced 

49 RNA binding affinity. In addition, we show that SAM promotes the formation 

50 of the active nsp10/nsp16 complex. Conversely, after cap methylation, the 

51 reaction turnover is speeded up by cap-1 RNA release and nsp10/nsp16 

52 complex dissociation, at the low intracellular SAH concentration. These 

53 results suggest that SAM/SAH balance is a regulator of the 2’-0- 

54 methyltransferase activity and raises the possibility that SAH hydrolase 

55 inhibitors might interfere with CoV replication cycle. The enzymatic and the 

56 RNA binding assays developed in this work were also used to identify 

57 nsp16 residues involved in cap-0 RNA recognition and to understand the 

58 action mode of known methyltransferase inhibitors. 
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59 Introduction: 

60 Middle East respiratory syndrome coronavirus (MERS-CoV) is an emerging 

61 zoonotic betacoronavirus that was initially isolated from a 60-year-old Saudi 

62 Arabian man in September 2012 (1). Since then, MERS-CoV is steadily 

63 spreading in the Arabian Peninsula and contained secondary outbreaks have 

64 occurred in Europe, Africa, Asia and North America with more than 1728 

65 confirmed infected patients up to April 2016 and a 36% fatality rate (2). 

66 MERS-CoV infection may be asymptomatic or result in clinical symptoms 

67 ranging from mild respiratory illness to severe acute pneumonia, renal 

68 failure, systemic disorder (1) and severe neurological syndrome (3). MERS- 

69 CoV belongs to the lineage C of the genus betacoronavirus and is 

70 phylogenetically related to the bat coronaviruses HKU4 and HKU5 (4-7). 

71 The bats are the host of a variety of coronaviruses including 

72 betacoronavirus genetically very closely related to the MERS-CoV which 

73 crossed the species barrier into dromedary camels (2). The virus was 

74 detected in Camelus dromedarius milk and nasal swabs, suggesting that 

75 camels are the probable source for zoonotic transmission of the virus to 

76 humans (6, 8). Recent works found at least five lineages of MERS-CoV in 

77 camels and identified six recombination events in MERS-CoV which may 

78 raise the virus’ pathogenicity (9). Human-to-human transmission requires 

79 close contact and occurred principally in health care settings (8)(10). 

80 Currently there are no approved antiviral treatments or vaccines available 

81 against MERS-CoV infection. 

82 Following binding of MERS-CoV spike envelope proteins to the dipeptidyl 

83 peptidase-4 receptor (DPP4, also known as CD26) (11), the virus genome 
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84 is released in the host cytoplasm. MERS-CoV genome is a polycistronic 

85 positive-stranded RNA of 30119 nt in length organized in eleven open 

86 reading frames (ORFs). The 3’ part of the genome contains at least nine 

87 ORFs that encode structural and virus accessory proteins, which are 

88 translated from a set of subgenomic RNAs. The two large ORFs (la and b) 

89 that encode non-structural proteins (nsp) are in the 5’ proximal two thirds of 

90 the genome (12). After viral infection, the 5’ region of capped genomic RNA 

91 is directly translated into the polyprotein ppla and pplab, which are the 

92 precursors of the CoV nsps. The pplab protein involves a ribosomal 

93 frameshift during RNA translation. These precursors proteins are cleaved 

94 by viral proteinases into 11 and 16 nsps, respectively, that form the 

95 replication-transcription complex (RTC) (12, 13). RTC harbours the catalytic 

96 activities required for the synthesis of genomic and subgenomic RNAs and 

97 contains most of the enzymes involved in the formation of the cap structure 

98 that decorates the 5’ end of viral mRNA. 

99 Capping is a common modification of the 5’ end of eukaryotic mRNA. The 

100 cap structure consists of a methylated guanosine at position 7 linked to the 

101 first transcribed nucleotide by a 5’-5’ triphosphate bridge ( 7m GpppN). The 

102 cap is co-transcriptionally added to the 5’ end of nascent mRNA after 

103 synthesis of the first 20-30 nucleotides (11,12). Cap structure synthesis 

104 involves four sequential reactions catalysed by an RNA 5’ triphosphatase 

105 (TPase), a guanylyltransferase (GTase), a guanine N7 methyltransferase 

106 (N7-MTase) and a 2’-0-MTase (15). The cap structure plays several key 

107 roles in mRNA functionality. First, capping is essential for mRNA stability 

108 and limits its degradation by cellular 5’-3’ exoribonucleases, such as XRN1 
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109 (16). In the cytoplasm, the cap-0 structure ( 7m GpppN) ensures efficient 

110 mRNA recognition by the eukaryotic translation Initiation Factor 4E (elF4E) 

111 (15). RNA capping also provides a molecular signature for the 

112 discrimination between self and non-self RNA. Indeed, viral double stranded 

113 RNA, 5’-triphosphate RNA and incorrectly capped RNA are detected by 

114 intracellular pathogen recognition receptors (PRRs). Among these PRRs, 

115 Retinoic acid-inducible Gene (RIG)-like receptors, such as RIG-1 and 

116 Melanoma differentiation-associated protein 5 (Mda5), can detect uncapped 

117 5’-triphosphate RNA and also cap-0 RNA (17-20). Detection of erroneously 

118 capped RNA in the cytoplasm induces a signal transduction cascade that 

119 initiates an antiviral response through interferon production (21). Among the 

120 interferon-stimulated genes (ISG), interferon-induced protein with 

121 tetratricopeptide repeats 1 (IFIT 1) can recognize mis-capped RNA and 

122 inhibit its translation (22, 23). 

123 Like other CoVs, MERS-CoV replicates in the cytoplasm of infected cells 

124 where it should produce its own capping machinery. CoVs seem to adopt 

125 the eukaryotic canonical capping pathway with four sequential events that 

126 involve several viral nsps: (i) first, the 5’-3’helicase/nucleoside 

127 triphosphatase (NTPase) nsp13 hydrolyses the y-phosphate from nascent 

128 5’-triphosphate RNA (24); (ii) then, a guanosine monophosphate (GMP) 

129 molecule is transferred to 5’-diphosphate RNA by a yet-unknown 

130 guanylyltransferase, thus forming a primitive cap structure (GpppN); (iii) 

131 then, the guanosine is methylated, at the N7 position, by nsp14 in the 

132 presence of the methyl donor S-adenosyl-L-methionine (SAM) to produce 

133 the cap-0 structure ( 7m GpppN) and S-adenosyl-homocysteine (SAFI) as by- 
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134 product (25); (iv) finally, the nsp10/nsp16 complex methylates the 2’-OH 

135 group (ribose) of the first transcribed nucleotide of the viral RNA, leading to 

136 the conversion of the cap-0 into a cap-1 structure ( 7 m GpppN 2 o m) (26, 27). 

137 Thus, in CoVs, RNA cap methylation might involve at least three proteins 

138 (nspIO, nsp14 and nsp16) and it is known that both nsp14 and nsp16 can 

139 interact with nspIO (28). 

140 CoV 2’-0-MTase activity is mediated by nsp16 that contains both a 

141 conserved K-D-K-E catalytic tetrad, which is characteristic of SAM- 

142 dependent 2’-0-MTases, and a conserved SAM-binding site (29). CoV 2’-0- 

143 MTase activity was first demonstrated using in vitro biochemical assays with 

144 purified feline CoV (FCoV) nsp16 (30). However, recombinant SARS-CoV 

145 nsp16 was inactive using similar experimental conditions. It was then 

146 demonstrated that SARS-CoV nsp16 needs to interact with nspIO to 

147 become catalytically active (26). The nsp10/nsp16 complex MTase activity 

148 was demonstrated using synthetic capped N7-methylated RNA and longer 

149 RNAs that mimic the 5’ end of the SARS-CoV genome (26). In contrast, 

150 RNA with an unmethylated cap structure (Gppp-RNA) was not recognized 

151 by the nsp10/nsp16 complex and no enzymatic activity was detected. 

152 Several mutagenesis studies of SARS-CoV nspIO and nsp16 confirmed the 

153 key role of the K-D-K-E catalytic tetrad for the 2’-0-MTase activity (31) and 

154 demonstrated that the interaction between nspIO and nsp16 is absolutely 

155 required for this activity (31, 32). The molecular basis governing nsp16 

156 stimulation by nspIO was partially elucidated by the crystal structure 

157 determination of the SARS-CoV nsp10/nsp16 complex (27, 31). Nsp16 

158 adopts the catechol O-MTase fold containing alternating p strands (P1-P7) 
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159 and a helices (aZ and aA- aE) that form a seven-stranded |3 sheet with 

160 three a helices on each side (31, 33). In addition, nspIO binds to nsp16 

161 through a 930 A 2 activation area on nspIO and stabilizes nsp16 (26, 27, 31, 

162 34). Moreover, structural and biochemical analyses also suggest that nspIO 

163 binding extends and narrows the RNA binding groove to accommodate the 

164 RNA substrate and enhances nsp16 RNA- and SAM-binding properties. 

165 Although cap 2’-0-MTases are present in different virus families (16), the 

166 exact role of this RNA-cap modification was only recently unravelled 

167 through reverse genetic studies performed, among others, on CoV (17-20). 

168 Single point mutations in the conserved K-D-K-E tetrad of nsp16 have 

169 revealed a barely attenuated phenotype in infected cells (17, 35). In 

170 contrast, infection of small animal models by viruses expressing nsp16 with 

171 active-site substitution showed a robust reduction of viral titers and less 

172 severe infection symptoms (weight loss, lung histology and breathing 

173 function) concomitantly with a strong antiviral response, possibly linked to 

174 stimulation of the innate immunity. Indeed, it was shown that incompletely 

175 capped RNA can be detected by RIG-1 or Mda5 (17, 19, 20), leading to the 

176 initiation of a signalling cascade that stimulates the secretion of type-1 IFN 

177 and induces an antiviral response in neighbouring cells. Altogether these 

178 observations indicate that the inhibition of the viral 2’-0-MTase activity 

179 might stimulate the detection of viral RNAs by RIG-1 or Mda5 and restore 

180 the host antiviral response mediated by ISG such as IFIT proteins (36). 

181 In this study, we assessed the biochemical activity of the MERS-CoV 

182 nsp10/nsp16 complex MTase. We show that MERS-CoV nsp16 possesses 

183 2’-0-MTase activity. As previously shown nsp16 is stimulated by nspIO, 
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184 and leads to 2’-0-methylation of cap-0 RNA (34). Biochemical assays were 

185 used to decipher the fine-tuning of the reaction turnover. Mutagenesis 

186 combined with RNA binding assays identified the residues essential for cap- 

187 0 recognition and those forming the RNA binding groove. Finally, using our 

188 radioactive enzymatic assay we screen a small subset of MTase inhibitors 

189 blocking nsp10/nsp16 MTase activity and the compounds blocking RNA 

190 recognition were identified using the fluorescence polarization assay. 
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191 Materials and methods 

192 Plasmid constructs 

193 The expression vectors pDEST-14/6 His-nspIO, pDEST-14/6 His-nsp16 

194 containing the coding sequence of the human betacoronavirus 2c 

195 EMC/2012 (GenBank: JX869059.2), a MERS-CoV strain, were provided by 

196 Eric Snijder’s team. For alanine scanning, nsp16 mutants were generated 

197 by single site amino acid substitution to alanine using the pDEST-14/6His- 

198 nsp16 plasmid and the Quickchange site directed mutagenesis kit (Agilent), 

199 according to the manufacturer’s instructions (mutagenic primers are listed in 

200 Table SI). The mutant numbering starts at the beginning of the nsp16 

201 sequence: Y30A, K31A, R38A, H41A, N43A, K46A, Y47A, K76A, D99A, 

202 N101A, D130A, Y132A, D133A, T136A, K137A, F149A, F150A, K170A, 

203 H174A, Y181A, E203A. All constructs were verified by DNA sequencing 

204 (Eurofins MWG operon). 

205 

206 Expression and purification of the MERS-CoV nspIO, nsp16 and 

207 nsplO/nspI6 proteins 

208 MERS-CoV nspIO and nsp16 fusion proteins (N-terminal hexahistidine tag) 

209 were expressed in E. coli C2566 cells that contain the pRARE-2 plasmid. 

210 Transformed bacteria cells were grown at 37°C in LB and 2YT medium, 

211 containing 100 pg/ml ampicillin and 17 pg/ml chloramphenicol. Protein 

212 expression was induced by addition of 0.5 mM isopropyl |3-D- 

213 thiogalactopyranoside (IPTG). After overnight incubation at 17°C, cells 

214 expressing nspIO, nsp16 or a mixture of nspIO and nsp16 (equal volume) 

215 were pelleted by centrifugation (13,000 x g, 10 min) and frozen before 
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216 resuspension in lysis buffer (50 mM Hepes pH 7.5, 300 mM NaCI, 30 mM 

217 imidazole, 10% glycerol supplemented with 1 mM PMSF, 0.25 mg/ml 

218 lysozyme and 10 pg/ml DNase I). After sonication and clarification (80,000 x 

219 g, 4°C, 30 min), the supernatant were incubated with HisPur™ Cobalt resin 

220 (Thermo Scientific) at 4°C with gentle shaking for 30 min. After washing in 

221 buffer W (50 mM Hepes pH 7.5, 40 mM imidazole, 10 % glycerol and 1 mM 

222 TCEP) containing 300 mM or 500 mM NaCI, bound proteins were eluted 

223 with buffer W supplemented with 250 mM imidazole. Recombinant nspIO 

224 and nsp16 proteins were then concentrated by ultrafiltration devices with 

225 molecular weight cut-offs (MWCO) of 5 kDa (Millipore) and 10 kDa 

226 (Sartorius), respectively. 

227 Nsp10/nsp16 complex was purified on a Superdex 200 column 16/60 (GE) 

228 equilibrated with buffer W containing 500 mM NaCI (GE AKTApurifier). After 

229 peak collection, fractions containing the nsp10/nsp16 complex were 

230 concentrated by ultrafiltration using 10 kDa Vivaspin 15 (Sartorius). All 

231 purified proteins were analysed by SDS-PAGE followed by Coomassie blue 

232 staining. Purified enzymes were stored at -20°C in 50% glycerol for 

233 enzymatic assays. 

234 

235 Methyltransferase assays 

236 Assay to test the MTase activity were carried out in a reaction mixture (40 

237 mM Tris-HCI, pH 8.0, 1 mM DTT, 1 mM MgCI 2 , 2 pM SAM and 0.33 pM 3 H- 

238 SAM (Perkin Elemer) in the presence of 0.7 pM synthetic RNA 

239 corresponding to the 5’ extremity of the MERS-CoV genome with various 5’ 

240 end modifications (tri-phosphorylated: pppGAUUUAA, cap: 
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241 GpppGAUUUAA, cap-0: 7m GpppGAUUUAA, or 2’-0-methylated cap: 

242 GpppG 20 mAUUUAA). Purified nspIO was added at a final concentration of 2 

243 pM in presence of 1.2 pM of nsp16 as the apparent interaction K d is about 

244 2 pM (Figure 2A). The purified nsp10/nsp16 complex was used at 1 pM. 

245 Reaction mixtures were incubated at 30°C and stopped at the indicated 

246 time points by diluting the reaction mixture with a ten-fold excess of 20 

247 microM ice-cold SAH (AdoHcy, New England Biolabs). Samples were then 

248 transferred to diethylaminoethylcellulose (DEAE) filters (Perkin Elmer) by 

249 using a Filtermat Harvester apparatus (Packard Instruments). 

250 Unincorporated 3 H-SAM was removed from the DEAE filters by several 

251 washes with 0.01 M ammonium formate, pH 8.0, H 2 0, and absolute ethanol. 

252 After drying, filters were incubated with BetaplateScint (Wallac) scintillation 

253 fluid before quantification of the 3 H methylation of the RNA substrates using 

254 a Wallac 1450 MicroBetaTriLux liquid scintillation counter (results were 

255 expressed as counts per minute, c.p.m.). 

256 For the MTase inhibition assays, 0.5 pM nsp10/nsp16 complex was mixed 

257 with 50 pM of each candidate inhibitor before addition of SAM (AdoMet, 

258 New England Biolabs) and 7m GpppGAUUUAA to start the reaction. The final 

259 DMSO concentration in the reaction mixtures was lower than 5%, and 

260 control reactions contained a similar DMSO concentration. Reaction 

261 mixtures were incubated at 30 °C for 30 min and then analysed by filter 

262 binding assay, as described above. 

263 The IC 50 values of SAH (AdoHcy), sinefungin, cap analogues ( 7m GpppA and 

264 7m GpppG) were determined with GraphPad Prism using a log (inhibitor) vs 

265 response variable slope equation. 
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266 

267 Fluorescence polarization assays 

268 RNA labelling: 10 pM RNA was labelled by ligation of 12.5 gM pCp- 

269 Cyanine-5 (pCp-Cy5) (Jena bioscience) at the 3’ extremity using 1 mM T4 

270 RNA ligase-1 (New England Biolabs) in 50 mM Tris-HCI pH 7.5, 10 mM 

271 MgCI 2 , 1 mM DTT and 1 mM ATP at 16°C overnight, according to the 

272 manufacturer’s instruction. T4 RNA ligase-1 was removed using Strata- 

273 Clean Resin (Agilent) and labelled RNA was collected after centrifugation in 

274 a Micro-Spin G-25 column (GE healthcare) to eliminate the free pCp-Cy5 

275 dye. 

276 Polarization assay: pCp-Cy5-labelled RNA was mixed with increasing 

277 concentrations of freshly purified nsp10/nsp16 complex, nspIO or nsp16 in 

278 binding buffer (50 mM Tris-HCI, pH 8, and 1 mM DTT) in a 384-well opaque 

279 microplate (Greiner bio-one). Fluorescence polarization (FP) was measured 

280 using a microplate reader (PHERAstar) with the optic module FP 

281 590/675 (590 nm excitation and 675 nm emission wavelengths). The 

282 dissociation constant values (K d ) were determined using non-linear 

283 regression analysis: site-specific binding with Hill slope equation (GraphPad 

284 Prism). The competitive effect of cap analogues on 7m GpppGAUUUAA-cy-5 

285 binding by 20 pM nsp10/nsp16 complex was determined using increasing 

286 concentrations of GpppG or 7m GpppG (New England Biolabs) in binding 

287 buffer. Data were plotted in GraphPad Prism 5.0 and a nonlinear sigmoidal 

288 dose response curve fitted in order to determine IC50 values of RNA 

289 binding inhibition. 

290 
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291 Interferometry for nsp10/nsp16 interaction 

292 NspIO biotinylation: nspIO was biotinylated (at room temperature for 30 

293 min) using EZ-Link NHS-PEG4-biotin (Thermofisher) at a molar ratio of 1:1 

294 in biotinylation buffer (50 mM Hepes pH 7.5, 150 mM NaCI, 10% glycerol, 1 

295 mM TCEP). Biotinylated nspIO was separated from free biotin using NAP-5 

296 columns equilibrated with biotinylation buffer. 

297 Octet analysis: Assays were performed using streptavidin sensors pre- 

298 incubated in 50 mM Hepes pH 7.5, 150 mM NaCI, 10 % glycerol, 1 mM 

299 TCEP and 0.5 mg/ml BSA. Biotinylated nspIO (100 nM) was loaded on 

300 streptavidin-coated biosensors for 5000s to reach 0.8 nm binding. Various 

301 nsp16 concentrations (0 to 7.8 pM) were used for the association step 

302 (300s) performed in 50 mM Hepes pH 7.5, 150 mM NaCI, 10% glycerol, 1 

303 mM TCEP and 0.5 mg/ml BSA. The dissociation step was followed for 500s. 

304 Real time nsp10-nsp16 interaction kinetics were recorded by bio-layer 

305 interferometry (BLI: Octet RED96). The K d , K on and K 0 ff values were 

306 calculated using the forteBio software after subtraction of the reference 

307 sensor from the sample traces. All curves were fitted using a 1:1 interaction 

308 model. The steady state curve of nsp10/nsp16 interaction was traced using 

309 one site-specific binding with the Hill slope equation (GraphPad Prism). 

310 Synthesis of RNA substrates 

311 RNAs were chemically synthesized on a solid support using an ABI 394 

312 synthesizer. After RNA elongation with 2’-0-pivaloyloxymethyl phosphoramidite 

313 monomers (Chemgenes, USA) (37), the 5’-hydroxyl group was phosphorylated 

314 and the resulting H-phosphonate derivative oxidized and activated into a 
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315 phosphoroimidazolidate derivative to react with pyrophosphate (pppRNA) (38) 

316 or guanosine diphosphate (GpppRNA) (39). After deprotection and release from 

317 the solid support, pppRNAs and GpppRNAs were purified by IEX-HPLC and 

318 their purity (>95%) confirmed by MALDI-TOF spectrometry. N7-methylation of 

319 purified GpppRNAs was then performed enzymatically using an N7-hMTase 

320 (39). 

321 

322 Mass spectrometry and IEX-HPLC analysis 

323 Crude pppRNAs and GpppRNAs were analysed and purified by IEX-HPLC on a 

324 Dionex Ultimate 3000 apparatus equipped with an UV detector at the 

325 wavelength of 260 nm and an anion exchange DionexDNAPac® PA200 column 

326 (4 x 250 mm for analysis or 9 x 250 mm for semi-preparative purpose) with a 

327 flow rate of 1.5 ml/min for analysis or 5 ml/min for semi-preparative purposes. 

328 The following eluents were used: Buffer A (20% CH 3 CN in 25 mM Tris-HCI, pH 

329 8), Buffer B (20% CH 3 CN containing 200 mM NaCI0 4 in 25 mM Tris-HCI, pH 8). 

330 After purification, pure fractions containing the desired RNA were pooled in a 

331 100 ml round-bottomed flask and concentrated to dryness under reduced 

332 pressure. Residues were desalted on a Cis cartridge (Sep-Pak® Classic). 

333 MALDI-TOF mass spectra were recorded on a Voyager-DE (PerSeptive 

334 Biosystems, USA) or on an AXIMA Assurance spectrometer (Shimadzu corp., 

335 Japan) equipped with a N 2 laser (337nm) using 2,4,6-trihydroxyacetophenone 

336 as saturated solution in a mixture of acetonitrile/0.1M ammonium citrate solution 

337 (1:1, v/v) for the matrix. The ionization mode was a linear negative ion mode. 

338 Analytical samples were mixed with the matrix in a 1:5 (v/v) ratio, crystallized on 

339 a 100-well stainless steel plate and analysed. 
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340 Results 

341 MERS-CoV nsp16 assembles with nspIO to constitute an active SAM- 

342 dependent cap-0 2’-0-methyitransferase. 

343 The CoV nsp16 2’-0-MTase contains both a conserved K-D-K-E catalytic 

344 tetrad, which is characteristic of SAM-dependent 2’-0-MTase, and a 

345 conserved SAM-binding site (29). While purified feline CoV (FCoV) nsp16 

346 shows a weak 2’-0-MTase activity in absence of nspIO, SARS-CoV nsp16 

347 2’-0-MTase activity is strictly dependent on its interaction with nspIO (30, 

348 34). MERS-CoV nsp16 shows -65% and -53% amino acid residue 

349 similarity with SARS-CoV and FCoV nsp16, respectively; and requires 

350 nspIO for nsp16 2’-0-MTase activity (34). To assess the mechanism of fine- 

351 tuning of MERS-CoV nsp16 2’-0-MTase activity by nspIO, recombinant 

352 nspIO and nsp16 with an N-terminal His 6 -tag were produced and purified, 

353 and SDS-PAGE analysis showed that nspIO and nsp16 migrating at the 

354 expected molecular weight (Figure 1A, 16 kDa and 37 kDa, respectively) 

355 and there identity was confirmed by mass spectrometry after trypsin 

356 digestion (not shown). The nsp10/nsp16 complex was obtained by co-lysis 

357 of an equal volume of His6-nsp10 and His6-nsp16 expressing bacteria and 

358 purification by metal affinity chromatography (IMAC) and size-exclusion 

359 chromatography. The chromatogram (Figure IB) indicated that the 

360 nsp10/nsp16 complex was detected as one major peak eluting at a volume 

361 corresponding to an nsp10/nsp16 heterodimer (~ 50 kDa). 

362 Then, in w'froMTase assays were performed using RNA oligonucleotides 

363 with sequences corresponding to the 5’-end of the MERS-CoV genome with 


D 

o 


Downloaded from http://jvi.asm.org/on January 2, 2017 by UCLA BIOMEDICAL LIB/SERIALS 


cfWology JVI Journal of Virology Accepted MdnUSCfTpf Posted OfTlme 


364 various 5’-end modifications (pppGAUUUAA, GpppGAUUUAA, 

365 7 m GpppGAUUUAA, GpppG 2 o mAUUUAA, 7 m GpppG 2 ' 0 mAUUUAA) in the 

366 presence of nspIO or nsp16 alone or as a complex and of 3 H-SAM as 

367 methyl donor. Quantification of the amount of [ 3 H]-CH 3 transferred to RNA 

368 over time using a DEAE filter-binding assay indicated the absence of 

369 MTase activity when nspIO or nsp16 alone was used (Figure 1C). 

370 Conversely, RNA carrying a cap-0 structure ( 7m GpppGAUUUAA) was 

371 methylated in a time-dependent manner when incubated with both nsp16 

372 and nspIO or with the purified nsp10/nsp16 complex. As MTase activity was 

373 detected only using 7m GpppGAUUUAA but not GpppG 2o mAUUUAA or 

374 7 m GpppG 2 omAUUUAA, we hypothesized that the MERS-CoV MTase 

375 converts cap-0 into cap-1 RNA. The addition of a single methyl group on 

376 7m GpppGAUUUAA was confirmed by comparing the molecular weight of 

377 7m GpppGAUUUAA before and after incubation with nsp10/nsp16 and SAM. 

378 We detected a mass increase of about =14.2 Da by MALDI-TOF mass 

379 spectrometry analysis indicating that one methylation event has occurred. 

380 We conclude that MERS-CoV nsp16 carries the cap 2’-0-MTase activity 

381 and is stimulated by nspIO. The nspIO stimulation effect on nsp16 was 

382 confirmed by assessing nsp16 2’-0-MTase activity in the presence of 

383 increasing concentrations of nspIO. An apparent dissociation constant for 

384 the nsp10/nsp16 complex could be derived from the stimulation curves (K d 

385 = 2 ± 0.1 pM, Figure 2A). Moreover, the nsp10/nsp16 MTase activity 

386 required the addition of divalent ions (Mg ++ or Mn ++ ) and was inhibited by 

387 EDTA (Figure 2B). The optimum of activity peaked between pH 8 and 8.5 

388 (Figure 2C). The optimal concentration of methyl-donor (SAM) (K d = 4.3 ± 
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389 0.6 |jM) was also determined by measuring the MTase activity of the 

390 complex in the presence of increasing SAM concentrations (Figure 2D). 

391 Finally, MTase assays performed using the optimal conditions deduced 

392 from these biochemical assays and cap-0 RNAs with sequences 

393 corresponding to the 5’-end of various viruses showed that the nsp10/nsp16 

394 complex methylated cap-0 RNAs that started with either an A or a G (Figure 

395 2E). Taken together, these results suggest that MERS-CoV nsp16 acts as 

396 an nspl 0-dependent 2’-0-MTase that converts cap-0 into cap-1 RNA 

397 structures with no obvious specific recognition mechanism of the MERS- 

398 CoV sequence. Moreover, the absence of 2’-0-MTase activity using 

399 GpppRNA suggests that nspl 6 can discriminate cap-0 structure from an 

400 unmethylated cap structure. 

401 

402 The MERS-CoV nsp10/nsp16 complex recognizes specifically a cap-0 

403 structure; and downstream RNA nucleotides promote loading of cap-0 

404 substrate and their 2’0 methylation 

405 RNA binding assays were then developed to determine why the 

406 nsp10/nsp16 complex was active only on cap-0 RNA. For this purpose, the 

407 3’-end of various 5’-end-modified RNA was labelled with pCp-Cy5 and then 

408 their binding to freshly purified MERS-CoV nsp10/nsp16 complex was 

409 measured in fluorescence polarization assays. The nsp10/nsp16 complex 

410 recognized only cap-0 RNA ( 7m GpppRNA) with a K d of -0.35 pM (Figure 3A). 

411 No significant binding was detected using the other RNAs (pppRNA, or 

412 GpppRNA, GpppG 2 ’omRNA, 7m GpppG 2o mRNA). This observation is 
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413 consistent with the MTase assay results (Figure 1C) and indicates that the 

414 absence of nsp10/nsp16 2’-0-MTase activity using GpppRNA is due to lack 

415 of recognition of such capped RNA. It also suggests that the nsp10/nsp16 

416 complex has an RNA binding site that recognizes specifically cap-0 RNA. 

417 Moreover, the ability of this binding site to recognize cap-0 RNA is strongly 

418 regulated by the 2’-0-methylation status of the substrate RNA. Indeed, 

419 7m GpppG 2 'omRNA, which corresponds to the nsp16 2’-0-MTase reaction 

420 product, did not bind to the nsp10/nsp16 complex (Figure 3A). This 

421 indicates that cap-1 RNA generated by the nsp10/nsp16 2’-0-MTase 

422 activity is actively released, promoting the reaction turnover. 

423 As nsp16 2’-0-MTase activity strictly depends on its interaction with nspIO 

424 (Figures 1C and 2A), we then asked whether nspIO regulates nsp16 RNA- 

425 binding properties. Comparison of nspIO, nsp16 and nsp10/nsp16 complex 

426 binding to cap-0 RNA ( 7m GpppRNA) using a fluorescence polarization assay 

427 showed that nspIO barely interacted with cap-0 RNA (Figure 3B). 

428 Conversely, both nsp16 and nsp10/nsp16 complex bound to cap-0 RNA 

429 with an apparent K d of ~1.4 and -0.33 pM, respectively. The polarization 

430 signal detected after nsp16 binding is higher than that obtained in the 

431 presence of the nsp10/nsp16 complex. This suggests that the RNA might 

432 bind to nsp16 dimers. Our result indicates that only nsp16 harbours the cap- 

433 0 binding site and that RNA binding is enhanced by nsp16 interaction with 

434 nspIO. As nsp16 alone binds cap-0 RNA (Figure 3B) but is not methylated 

435 (Figure 1C), we suggest that nspIO is an allosteric regulator of nsp16 2’-0- 

436 MTase activity. In addition, the methyl donor SAM and, to a lesser extent, 

437 the reaction by-product SAH increased the affinity for the RNA substrate 
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438 (Figure 3C). Conversely, magnesium ions (Mg 2+ ), which stimulate 

439 nsp10/nsp16 MTase activity (Figure 2B), did not seem to affect RNA 

440 binding (Figure 2F) in these experimental conditions, indicating that 

441 magnesium ions are not required for RNA recognition. 

442 The MTase assay showed that the nsp10/nsp16 complex is able to 

443 methylate cap-0 RNA of variable RNA sequences (Figure 2E), starting with 

444 either m GpppA or m GpppG. Flowever, the MTase activity is weaker using 

445 short RNA substrates (Figure 3E). Therefore, we tested whether 

446 nucleotides downstream the cap structure contributed to nsp10/nsp16 

447 complex binding to RNA. To this aim, fluorescence polarization assays were 

448 performed using 3’-end labelled RNA sequences of increasing length that 

449 mimicked the 5’-end of the MERS-CoV genome sequence. Quantification of 

450 their interaction with the nsp10/nsp16 complex (Figure 3D and Table 1) 

451 showed that the cap analogue 7m GpppG-cy5, which was not methylated by 

452 nsp10/nsp16 (Figure 3E), barely interacted with the nsp10/nsp16 complex. 

453 In contrast, a significant interaction was observed using 7m GpppGAUU-cy5. 

454 The affinity increased then with the RNA length and the optimal affinity was 

455 observed with 7m GpppGAUUUAAGUG-cy5 (Figure 3D). The MTase activity 

456 measured using a filter-binding assay followed the same trend (Figure 3E). 

457 Identification of nsp16 residues playing a key role in RNA recognition 

458 and MTase activity 

459 The crystal structure of the SARS-CoV nsp10/nsp16 complex with SAH has 

460 been reported at 2.0 A resolution (31). As no structure with a cap analogue 

461 or substrate RNA is available yet, the nsp16 residues involved in RNA 

462 binding were tentatively inferred by modelling the RNA in the catalytic site 
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463 (Figure 4: structural model based on the SARS-CoV and vaccinia MTase 

464 VP39 structures). Based on the structural model, alanine mutagenesis was 

465 performed and the effects of nsp16 mutations on MTase activity (Figure 5B) 

466 and on RNA binding properties (Figure 5A) were determined by filter 

467 binding assay and fluorescent polarization assay, respectively. Alanine 

468 substitution of the conserved residues in the catalytic tetrad (K-D-K-E) 

469 almost completely abolished the 2’-0-MTase activity, as expected. With the 

470 exception of D130A, these mutations also strongly reduced RNA binding, 

471 indicating that these catalytic residues also participated in RNA binding 

472 process. Mutation of the residues adjacent to the putative SAM binding site 

473 (H41, N43, D99, F149 and F150) drastically reduced MTase activity, but 

474 had no significant effect on RNA binding. Mutation of K46, K170, E203, 

475 D133, R38, Y47 or Y181 concomitantly reduced MTase activity (80-100%) 

476 and RNA binding. These residues are localized in the putative RNA binding 

477 groove of nsp16 (Figure 4). In addition, mutation into alanine of Y30 and 

478 Y132, which are within two mobile a helices (26-38) and (130-148) in close 

479 proximity with the cap in the nsp10/nsp16 model (Figure 4), and of HI74, 

480 which is close to the 130-148 loop, also decreased nsp16 RNA binding 

481 properties and MTase activity. These observations confirm that the cap-0 

482 RNA binding site involves these three aromatic residues. 

483 SAM and SAH levels regulate nsp10-nsp16 interaction 

484 Nsp16 2’-0-MTase activity requires interaction with nsplO. To determine 

485 how co-substrate (SAM) availability might regulate the RNA-capping 

486 reaction pathway, the nsp16 interaction with nsplO was analysed using 
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487 biolayer interferometry. NspIO was biotinylated before immobilization on a 

488 streptavidin-coated biosensor. Binding of nsp16 was then monitored in the 

489 presence of various SAM/SAH concentrations, and the apparent affinity 

490 constants were determined (Figure 6A and Table 2). Nsp16 interaction with 

491 biotinylated nspIO increased in the presence of 100 pM SAM or SAH. The 

492 interaction kinetics was characterized by the same association rate (same 

493 k on ) in the presence or in the absence of SAM or SAH. Conversely, the 

494 addition of SAM or SAH increased nsp10/nsp16 complex stability by 

495 decreasing the dissociation kinetics (k 0 ff). This suggests that both SAM and 

496 SAH induce nsp10/nsp16 complex assembly. To further analyse how SAM 

497 or SAH regulated nsp10-nsp16 interaction, the dissociation of the 

498 nsp10/nsp16 complex, formed in the presence of 100 pM SAM, was 

499 assessed using dissociation buffers containing various concentrations of 

500 SAM or SAH (Figure 6B). The dissociation rate of the nsp10/nsp16 complex 

501 depended on SAM or SAH concentration. In the presence of 100 pM SAM, 

502 a concentration that mimics its intracellular level (40), the nsp10/nsp16 

503 complex stability was highest, as indicated by the slow nsp16 release. In 

504 contrast, after MTase reaction, which converts SAM into SAH, the latter is 

505 released from the enzyme at the lower intracellular concentration of SAH 

506 (~20 pM (41)). In turn, the release of SAH will speed-up the dissociation of 

507 the nsp10/nsp16 complex. Altogether, these data indicate that nspIO and 

508 nsp16 form a dynamic complex the association of which is promoted by 

509 SAM binding. The complex will dissociate when SAH, the reaction by- 

510 product, is released, thus favouring the reaction turnover. 

511 
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512 MERS-Co V nsp 16 MTase activity inhibition. 

513 It was previously shown that 2’-0-methylation of viral RNA cap structures 

514 limited the infected host’s antiviral response (42), suggesting that molecules 

515 blocking nsp16 2’-0-MTase activity might favour viral clearance in infected 

516 animals (36). In MTase inhibition assays, we tested 18 molecules known to 

517 block MTase activity (26, 43). Four of these compounds (Figure 7A) 

518 inhibited nsp10/nsp16 2’-0-MTase activity (between 45 and 100%) at a 

519 concentration of 50 pM. Among them, two SAM analogues (sinefungin and 

520 SAH: number 2 and 3 in Figure 7A) and showed the most potent inhibition. 

521 Cap-0 analogues ( 7m GpppG and 7m GpppA: number 18 and 16) also 

522 inhibited nsp16 2’-0-MTase activity (40-45%), whereas GpppG and GpppA 

523 (number 17 and 15) did not. This suggests that cap analogues block the 

524 binding of the RNA substrate in nsp16 RNA or cap binding site, in 

525 agreement with the weak interaction (Kd >20 pM) of these compounds with 

526 the nsp10/nsp16 complex (Figure 3D). 

527 To better define the 2’-0-MTase activity inhibition of sinefungin, SAFi, 

528 7m GpppG and 7m GpppA, the nsp10/nsp16 complex was pre-incubated with 

529 increasing concentrations of each inhibitor, and then, the MTase reaction 

530 was started by addition of 3 FI-SAM and the RNA substrates. The dose- 

531 response curves (Figure 7B) indicated that sinefungin, and SAH inhibited 

532 nsp10/nsp16-mediated 2’-0-MTase activity with IC 50 values in the 

533 micromolar range (7.4 and 7.0 pM, respectively, Table 3). The IC 50 values 

534 of the cap analogues 7m GpppG and GpppG (45 and 274 pM, respectively) 

535 (Table 3) suggest that 7m GpppG is a competitive inhibitor that blocks 

536 specifically the RNA cap-0 binding site. This was confirmed by measuring 
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537 the effect of these compounds on nsp10/nsp16 RNA binding activity. 

538 7m GpppG blocked RNA binding with an IC 50 value of 45 pM, which is in the 

539 same range of that of the inhibition of the MTase activity (Figure 7C). 

540 Conversely, the SAM analogues (SAH and sinefungin), which are supposed 

541 to enter in the nsp16 SAM binding site, barely affected RNA binding 

542 properties, as expected (Figure 7D and Table 4). Thus, the radioactive 

543 assay provides an efficient method to screen molecules to identify inhibitors 

544 of SAM-dependent MTase activity of the MERS-CoV nsp10/nsp16 complex. 

545 The fluorescence polarization assay developed here represents a tool to 

546 analyse their mode of action. 

547 
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548 Discussion: 

549 

550 MERS-CoV is the most recently discovered zoonotic virus that causes 

551 pathology in humans. Upon cell infection, MERS-CoV starts its replication 

552 by translation of the ORF la and 1b to produce 16 nsps of the 

553 replication/transcription complex (RTC)(12). In this study, we characterized 

554 the MERS-CoV 2’-0-MTase activity of nsp16 and show that nspIO is a 

555 cofactor required for nsp16 2’-0-MTase enzymatic activity, as previously 

556 reported for SARS and MERS-CoV (26, 34). The substrate of the reaction is 

557 RNA with a cap-0 structure. This observation is consistent with RNA binding 

558 assays showing that N7-methylated cap RNA is recognized by the 

559 nsp10/nsp16 complex, whereas unmethylated RNA is not. Thus, MERS- 

560 CoV cap methylation follows an obligatory sequence in which 2’-0- 

561 methylation might occur after N7 methylation mediated by nsp14. 

562 Consequently, it is likely that RNA cap synthesis by MERS-CoV follows the 

563 canonical capping pathway observed in eukaryotic cells and in some 

564 viruses, such as SARS-CoV, dengue virus, and West Nile virus (21). 

565 Although the cap-0 structure gives RNA binding specificity, cap analogues 

566 barely bind to the nsp10/nsp16 complex at the micro molar concentration of 

567 7m GpppG-cy5 used in our assay, and are not a good substrate for the 2’-0- 

568 MTase. The presence of additional nucleotides downstream of the cap 

569 structure increases the substrate recognition and 2’-0-MTase activity. 

570 These nucleotides are thus required for substrate stabilization in the RNA 

571 binding groove. The optimal recognition is reached with ten nucleotides 

572 capped RNA. This is consistent with the model proposed for the SARS-CoV 
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573 nsp10/nsp16 complex with 7m GpppA-RNA in which the first nucleotides are 

574 held by nsp16 and the following nucleotides might be positioned in a RNA 

575 binding groove stabilized by nspIO (27, 31). In addition, it is likely that the 

576 MERS-CoV 2’-0-MTase recognizes and methylates cap-0 RNA with no 

577 obvious sequence specificity. We also observed that the nsp10/nsp16 

578 complex methylated similarly cap-0 RNA with A or G as the first nucleotide. 

579 This observation is in contrast to SARS-CoV 2’-0-MTase that preferentially 

580 binds and methylates RNA with A at Ni position (( 7m GpppA-capped RNA, 

581 (27)). Using synthetic RNAs that mimic the 5’-end of the MERS-CoV 

582 genome, we did not observe any additional methylation (N2) or internal 

583 methylation, as indicated by the absence of enzymatic activity on cap-1 

584 RNA or pppRNA sequences. The absence of internal methylation was also 

585 confirmed using longer (about 500 nucleotides) pppRNA substrates (not 

586 shown). This observation contrasts with findings on other viral MTases 

587 (flavivirus NS5-MTase) that can perform internal methylation of adenosines 

588 (44). Nevertheless, the absence of internal methylation with the 

589 nsp10/nsp16 complex is consistent with the fact that cap-1 RNAs or 

590 pppRNAs do not interact with the nsp10/nsp16 complex. 

591 

592 To further understand the mechanism whereby nspIO promotes the 

593 stimulation of nsp16 2’-0-MTase activity by nspIO, we determined whether 

594 nspIO is required for substrate recognition process. Our findings indicate 

595 that nsp16 can recognize cap-0 RNA even in the absence of nspIO as 

596 already reported for FCoV nsp16 (30). Nevertheless, nspIO increases 

597 nsp16 RNA binding properties of nsp16 (3- fold) and in turn the 2’-0-MTase 
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598 activity can be detected only in the presence of both proteins under this 

599 condition. These observations suggest that nspIO acts as an allosteric 

600 regulator of nsp16 2’-0-MTase activity rather than as an RNA binding 

601 module. A similar allosteric effect was reported for the vaccinia virus and 

602 human RNA N7-MTases. Indeed, in the vaccinia virus family, the D1 

603 catalytic subunit alone is unstable and inactive and its stability and MTase 

604 activity are enhanced by the presence of D12 through the increase of GTP, 

605 SAM and GpppA binding affinity (45, 46). Similarly, the human RNA N7 

606 MTase (RNMT) is allosterically activated by RNMT Activating Miniprotein 

607 (RAM) that stabilizes its structure and favours the recruitment of methyl 

608 donors (47). 

609 

610 As nspIO is not the main player involved in RNA binding, we mapped the 

611 nsp16 residues involved in cap-0 RNA recognition. We could not determine 

612 the structure of nsp10/nsp16 with RNAs or cap analogues, although we 

613 performed many attempts (not shown). Therefore, we used alanine 

614 scanning mutagenesis to identify the residues involved in cap-0 RNA 

615 specific recognition and nsp16 MTase activity. Our results suggest that the 

616 cap-0 RNA structure is stacked between two flexible loops (26-38) and 

617 (130-148) by the aromatic residues Y30, Y132 respectively and the H174 

618 residue in close proximity. The alanine scanning results also indicate that 

619 the RNA chain of the substrate might be held by K46, K170, E203, D133, 

620 R38, Y47 and Y181 because both RNA binding and MTase activity 

621 decreased upon mutation of these residues. Interestingly, three residues 

622 belonging to the catalytic tetrad (K46, K170 and E203) might directly 
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623 participate in RNA recognition. These results corroborate the model of 

624 MERS-CoV nsp10/nsp16 (Figure 4) in which the RNA was positioned using 

625 vaccinia MTase VP39 structures (PDB entry: 1AV6) (27, 31, 48). The 

626 mutagenesis study also allowed identifying residues that strongly reduce 

627 MTase activity with a weak effect on RNA recognition (D99, F149, F150, 

628 H41 and N43). On the basis of their position in the structural model, we 

629 suggest that these residues might participate in SAM binding. 

630 

631 We then demonstrated that SAM/SAFI balance plays an important role in 

632 nsp10-nsp16 interaction. The methyl donor SAM, strengthens the 

633 interaction between nspIO and nsp16 at the physiological intracellular 

634 concentration ((100 pM, (40)) and in turn, the RNA binding properties of 

635 nsp16 are enhanced. This suggests that SAM stabilizes or induces small 

636 conformational changes of the enzyme leading to an increase in both RNA 

637 affinity and methylation by nsp10/nsp16. This observation is also consistent 

638 with thermal shift assay indicating that the nsp10/nsp16 complex is 

639 stabilized in the presence of SAM (AT°m = 1 deg). During RNA methylation, 

640 SAM is converted into SAH. The SAH will be next released by nsp16 

641 stimulating the dissociation of the nsp10/nsp16 complex. In addition, the 

642 reaction product (cap-1 RNA) barely interacts with the nsp10/nsp16 

643 complex and is released. Overall, these observations suggest that cap-0 

644 RNA methylation is linked to the association and dissociation of the 

645 nsp10/nsp16 complex (Figure 8). These data also suggest that the 

646 SAM/SAFI balance in infected cells is a key factor for cap RNA methylation. 

647 As, SAH hydrolase is a key enzyme for the maintaining of low levels of SAH 
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648 in cells, inhibitors of cellular SAH hydrolases, might have an antiviral effect 

649 on CoV replication. Interestingly, the SAH hydrolase inhibitor 3- 

650 deazaneplanocin A, has previously been shown to inhibit several viruses 

651 (49). It was recently demonstrated to limit Ebola virus replication and to 

652 lead to IFN- a production in virus-infected mice (50, 51). The SAH 

653 hydrolase inhibitors may thus confer a selective antiviral activity that differ 

654 from one virus to another depending on their RNA methylation levels (52). 

655 The inhibition may result to a direct effect of SAH concentration increase on 

656 the viral MTase or to indirect effects on cellular mRNA cap methylation 

657 which may then be recognized as "non-self" and, thus trigger an IFN 

658 response. 

659 

660 The viral 2’-0-MTase was highlighted as potential antiviral target (42). 

661 Indeed, such RNA methylation could be a marker of "self”, thus avoiding 

662 detection of viral RNA by Mda5/RIG-I sensors (17, 19, 20). In addition, one 

663 interferon-stimulated gene products, IFIT1, sequesters 2’0-unmethylated 

664 viral RNA and thereby blocks their translation, which results in inhibition of 

665 viral replication (53). Thus, compounds blocking the viral 2’-0-MTase could 

666 limit viral replication in infected animals and elicit a strong antiviral 

667 response, favouring viral clearance. Here, we show that SAM analogues, 

668 such as SAH and sinefungin, limit 2’-0-MTase activity with IC 50 values in 

669 the micromolar range (7.0 and 7.4 pM, respectively) with no inhibition of 

670 RNA binding, as expected. Conversely, a N7-methylated cap analogue 

671 ( 7 m GpppG) inhibits both 2’-0-MTase activity and RNA recognition with an 

672 IC 50 value of 45 pM. This indicates that 7m GpppG acts as a competitive 
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673 inhibitor that blocks cap-0 RNA recognition by the nsp10/nsp16 complex. 

674 This was confirmed by the absence of inhibition in the presence of 

675 unmethylated cap analogues. Altogether, these results might help 

676 developing an efficient assay to screen compounds that inhibit the SAM- 

677 dependent MERS-CoV nsp10/nsp16 complex. The biochemical tools and 

678 the RNA binding assay allow determining whether inhibitors are competing 

679 directly with the RNA substrate. It remains to be tested whether targeting 

680 the 2’-0-MTase activity can limit viral replication in cells or animal models 

681 that harbour the Mda5/RIG-I antiviral pathway (34). 

682 
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877 Figures 

878 

879 Figure 1: MERS-CoV nsp16 2-O-MTase activity is promoted by nspIO 

880 Recombinant MERS-CoV nspIO, nsp16 and the nsp10/nsp16 complex were 

881 expressed in E. coti and purified by affinity chromatography on IMAC columns. 

882 A) After SDS-PAGE separation of purified proteins (0.2 pg), gels were stained 

883 with Coomassie blue. MW, molecular size markers. B) The nsp10/nsp16 

884 complex was purified by gel filtration using a S200-16/60 column. The elution 

885 chromatogram (OD: 280 nm) shows one main peak that eluted at 90 ml and 

886 corresponded to the nsp10/nsp16 complex. C) The methyltransferase activity 

887 (MTase) of nsp16 (1.2 pM) and nspIO (2 pM) alone or together and of the 

888 nsp10/nsp16 complex (1 pM) was determined by monitoring the transfer of 3 H- 

889 CH 3 from SAM to RNA oligonucleotides with sequences corresponding to the 

890 5’-end of the MERS-CoV genome with various cap modifications 

891 (pppGAUUUAA, GpppGAUUUAA, 7m GpppGAUUUAA GpppG 2 . om AUUUAA and 

892 7m GpppG 2 'omAUUUAA). Assays were stopped after 5, 30, 60 and 180 min and 

893 the radioactivity associated with the RNA was determined by DEAE-filter 

894 binding assay. The bar graph presents the mean value and the standard 

895 deviation of three independent experiments. 

896 

897 Figure 2: Biochemical characterization of the MERS-Co V nsp 10/nsp 16 2’- 

898 O-MTase 

899 

900 The 2’-0-MTase activity of MERS-CoV nsp16 was characterized by filter 

901 binding assay as described in Figure 1C. A) Nsp16 2’-0-MTase activity is 
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902 promoted by nsplO. Nsp16 (2 pM) was incubated with 7m GpppGAUUUAA (0.7 

903 pM) in the presence of increasing concentrations of nsplO. The apparent K d of 

904 nsp10/nsp16 interaction (2 ± 0.1 pM; mean value and standard error of the 

905 mean, SEM, of three independent experiments) was deduced by curve fitting 

906 using Hill slope equations. B) Effect of divalent ions on nsp16 2’-0-MTase 

907 activity. MTase assays were performed in the presence or not of divalent ions 

908 (MgCI 2 ). Inhibition of the MTase reaction by EDTA was assessed and the 

909 activity was recovered by addition of MgCI 2 or MnCI 2 . The methyl transfer was 

910 measured after 5, 15, 30 and 60 min at 30°C (mean value ± SD). C) Optimum 

911 pH of the 2’-0-MTase activity mediated by nsp10/nsp16. D). Nsp10/nsp16 

912 MTase activity in the presence of increasing concentrations of SAM. The 

913 apparent K d of SAM (4.3 ± 0.6 pM) was deduced by curve fitting using a Hill 

914 equation and the GraphPad Prism program. E) Measurement of nsp10/nsp16 

915 MTase activity on cap-0 RNA sequences that correspond to the 5’-end of 

916 various viruses: MERS-CoV ( 7m GpppGAUUUAAGUGAAUA), West Nile virus 

917 ( 7m GpppAGUAGUUCGCCUG), Dengue virus ( 7m GpppAGUUGUUAGUCUA), 

918 Ebola virus ( 7m GpppGAUGAAGAUUAAG) after 30 min incubation at 30°C (n=3, 

919 mean value ± SD) .F) Magnesium ions do not influence nsp10/nsp16 interaction 

920 with cap-0 RNA. 7m GpppGAUUUAA-cy-5 was incubated with increasing 

921 concentrations of the nsp10/nsp16 complex with or without 1 mM MgCI 2 . The 

922 interaction with the cy-5-labelled RNA was followed by fluorescent polarization, , 

923 and the apparent calculated K d values are indicated on the top of the graph. 

924 

925 Figure 3: RNA binding properties of MERS-CoV nsplO, nsp16 and the 

926 nsp10/nsp16 complex 
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927 A) Nsp10/nsp16 recognizes cap-0 RNA. Different short RNA oligonucleotides 

928 (pppGAUUUAA, GpppGAUUUAA, 7m GpppGAUUUAA, GpppG 2o mAUUUAA or 

929 7 m GpppG 2 omAUUUAA) corresponding to the 5’-end of the MERS-CoV genome 

930 were 3’-end-labelled with pCp-cyanine-5 and incubated with increasing 

931 concentrations of nsp10/nsp16 complex. The interaction between the 

932 nsp10/nsp16 complex and each RNA was followed by measuring the 

933 fluorescence polarization signal at 675 nm. The apparent affinity constant (K d ) 

934 was calculated by non-linear regression analysis with a Hill slope equation and 

935 is indicated on the top of the graph. B) Nsp16 and the nsp10/nsp16 complex 

936 recognize cap-0 RNA. 7 m GpppGAUUUAA-cy-5 was incubated with increasing 

937 concentrations of nspIO, nsp16 or nsp10/nsp16 complex. The interaction with 

938 cy-5-labelled RNA was followed by measuring the fluorescence polarization, like 

939 in panel “A”, and the calculated K d is indicated on the top of the graph. C) SAM 

940 and SAH enhance the interaction of the nsp10/nsp16 complex with cap-0 RNA. 

941 7 m GpppGAUUUAA-cy-5 was incubated with increasing concentrations of 

942 nsp10/nsp16 complex in the presence or not (buffer alone) of 100 pM SAM or 

943 20 pM SAH. The interaction with cy-5-labelled RNAs was followed by 

944 measuring the fluorescence polarization and the calculated K d is indicated on 

945 the top of the graph. D) RNA length effect on nsp10/nsp16 binding was 

946 analysed using viral native 7m Gppp RNAs of different lengths ( 7m GpppG, 

947 7m GpppGAUU, 7m GpppGAUUUAA, 7m GpppGAUUUAAGUG and 

948 7 m GpppGAUUUAAGUGAAUA) in the presence of increasing concentrations of 

949 nsp10/nsp16. The interaction with cy-5-labelled RNAs was followed by 

950 fluorescence polarization and the calculated K d values are indicated in Table 1. 

951 E) Kinetics of MTase activity of the nsp10/nsp16 complex (0.25 pM) on 
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952 7m GpppG-RNAs (0.7 pM) of increasing length was measured by filter binding 

953 assay after 2, 4, 8, 16, 32 and 64 min incubation at 30°C in the presence of 3 H- 

954 SAM. The curve presents the mean value and the standard deviation of three 

955 independent experiments. All the fluorescent polarization experiments were 

956 performed twice independently and the K d values of each condition were 

957 calculated using GraphPad prism (n=2, mean value ± SEM) using one site- 

958 specific binding equation with Hill slope. 

959 

960 Figure 4: 3D structural model of the nsp10/nsp16 complex interaction with 

961 7m Gppp-RNA 

962 The nsp10/nsp16 model was built by alignment of SARS-CoV (PDB entry 3R24) 

963 and MERS-CoV nsp10/nsp16 sequences using the Swiss model and the 

964 PyMOL software. The cap-0 RNA/nsp16 model was built by aligning the MERS- 

965 CoV nsp16 sequence with that of the vaccinia MTase VP39 (PDB entry 1AV6) 

966 using the PyMol software. The surface of nsp16 residues is in blue and that of 

967 nspIO in pink. The nsp16 catalytic residue D130 is depicted in orange, while the 

968 catalytic and putative RNA binding residues K46, K170, E203 are in red. The 

969 putative cap-0 RNA binding site (light green) is delimited by Y30 and Y132, 

970 which are localized in the two mobile a helices (26-38) and (130-148) 

971 respectively, and by HI74 close to the last helix. The residues of the putative 

972 RNA binding groove are in cyan (Y181, D133, Y47 and R38). The SAM binding 

973 pocket involves the residues N43, D99, F149, F150 and H41 (in magenta); SAM 

974 is in grey and cap-0 RNA in yellow. The zoom focuses on the RNA binding 

975 groove and the SAM binding domain (ribbon structure). 

976 
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977 Figure 5: MERS-CoV nsp16 mutations that affect RNA recognition and/or 

978 MTase activity 

979 A) Alanine scanning mutagenesis. Mutations were introduced in the nsp16 

980 clone and after co-lysis, the nsp10/nsp16 complex was purified by affinity 

981 chromatography on IMAC. The affinity (K a = 1/kci) of each mutant protein for 

982 7m GpppGAUUUAA-cy-5 was measured by fluorescence polarization as 

983 described in Figure 3 The affinity (K a = 1/ K d) of each mutant protein for 

984 7m GpppGAUUUAA-cy-5 was measured by fluorescence polarization as 

985 described in Figure 2. The K a values were calculated from the Hill equation 

986 using GraphPad prism (n=2, mean value ± SD). The K a values were only 

987 estimated for the mutants indicated by an asterisk because the RNA binding 

988 plateau was not reached in the fluorescence polarization experiment. B) The 

989 MTase activity of each mutant (1 pM) was determined after incubation with 

990 7m GpppGAUUUAA at 30°C in the presence of 3 H-SAM for 60 min. The bar 

991 graph presents the mean value and the standard deviation of three independent 

992 experiments. 

993 

994 Figure 6: SAM and SAH modulate nsp10-nsp16 interaction 

995 A) Nsp10/nsp16 steady state assembly: real-time nsp16 binding to biotinylated 

996 nspIO was measured by Octet biolayer interferometry (BLI). Streptavidin 

997 biosensors coated with 100 nM biotinylated nspIO were used to measure nsp16 

998 (0 to 7.8 pM) association and dissociation. The reference sensor (0 pM nsp16) 

999 values were subtracted from the sample traces. All curves were fitted with a 1:1 

1000 model by using the Octet biolayer program and then the steady state curve was 

1001 traced using the site-specific binding equation on GraphPad Prism. The 
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1002 different K d values of nsp10-nsp16 interaction in the presence or not of 100 pM 

1003 SAM or 20 pM or 100 pM SAH are presented in Table 2. B) Effects of SAM and 

1004 SAH on the dissociation kinetics of nsp16 from immobilized nsplO. 7.8 pM 

1005 nsp16 was bound to sensor-immobilized biotinylated nsplO in the presence of 

1006 100 pM SAM. Then, the sensor was moved to dissociation buffer alone or 

1007 containing 5, 20 or 100 pM of SAM or SAH. Nsp16 binding to immobilized 

1008 nsplO was normalized for each trace. 

1009 

1010 Figure 7; Inhibition of nsp10/nsp16 complex MTase activity and RNA 

1011 binding 

1012 A) Bar graph showing the MTase inhibition activity (percentage) of each 

1013 candidate inhibitor (final concentration of 50 pM). The MTase assay was 

1014 performed as in Figure 1C by incubating 0.5 pM of nsp10/nsp16 complex with 

1015 0.7 pM 7m GpppGAUUUAA and 3 H-SAM at 30°C for 30min. The methyl transfer 

1016 to RNA was measured by filter binding assay in the absence (5% DMSO) or in 

1017 the presence of each inhibitor (n=2, mean value ± SD) 1: DMSO, 2: sinefungin, 

1018 3: SAH, 4: SIBA, 5: ribavirin, 6: ribavirin-TP, 7: N-(5-chloro-2-methoxyphenyl)- 

1019 3,4-dimethoxy-N-[2-(4-methyl-1-piperazinyl)-2-oxoethyl]benzenesulfonamide, 8: 

1020 2-[[4-benzyl-5-(3-pyridinyl)-4H-1,2,4-triazol-3-yl]thiol]-N-(5-methyl-1,3,4- 

1021 thiadiazol-2-yl) acetamide, 9: N-1-cyclopropyl-N-2-(2,5-dimethoxyphenyl)-N-2- 

1022 (methylsulfonyl) glycinamide, 10: 4-bromophenyl 3-(3,5-dioxo-4-azatetracyclo 

1023 [5.3.2.0-2.6-.0-8,10-] dodec-11-en-4-yl) propanoate, 11: 1'-[(4-tert- 

1024 butylphenyl)sulfonyl]-1,4'-bipiperidine-4'-carboamide, 12: N'-(3,4- 

1025 dichlorobenzylidene)-2-(2-pyridinylthio) acetohydrazide, 13: GTP, 14: 7m GTP, 

1026 15: GpppA, 16: 7m GpppA, 17: GpppG, 18: 7m GpppG. B) Dose response of the 
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1027 inhibitory effect of sinefungin, SAH and cap analogues on nsp10/nsp16 MTase 

1028 activity (MTase assays performed as in panel A; n=2, mean value ± SEM). The 

1029 IC 50 values deduced using GraphPad Prism and the log (inhibitor) vs response 

1030 variable slope equation are shown in Table 3. C) The competition effect of cap 

1031 analogues on 7m GpppGAUUUAA-cy-5 binding to the nsp10/nsp16 complex was 

1032 analysed by fluorescence polarization, as described in Figure 3. Measurements 

1033 were performed in the presence of increasing concentrations of cap analogues. 

1034 D) The effect of 20 pM SAH, 20 pM sinefungin or buffer alone on the 

1035 nsp10/nsp16 complex interaction with 7m GpppGAUUUAA-cy-5 was assessed by 

1036 fluorescence polarization, as described in Figure 3C. The K d values deduced 

1037 from the curve fitting are presented in Table 4 (n=2, mean value ± SEM). 

1038 

1039 Figure 8: Model of the nsp10/nsp16- 7m GpppRNA reaction turnover 

1040 regulated by the SAM/SAH balance 

1041 The nsp10/nsp16 model was built as described previously in Figure 4. The 3D 

1042 representation of nsp16 structure is in blue, 7m GpppGAAAAA in yellow and 

1043 nspIO in pink. According to our model, first, nsp16 binds to cap-0 RNA and 

1044 loads a SAM molecule (grey, 100 pM in infected cells), thus allowing the 

1045 recruitment of the nspIO allosteric activator. Then, nsp16 catalyses cap RNA 

1046 methylation and converts cap-0 into cap-1 RNA. Methylated RNA is released 

1047 from the nsp16 cap-binding site. As SAH, in orange, is supposed to be at lower 

1048 concentration than SAM in infected cells (20 pM vs 100 pM), SAH should be 

1049 released from the nsp10/nsp16 complex and favour its dissociation and the 

1050 reaction turnover. 

1051 
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1052 Table 1: K d values of the interaction of the nsp10/nsp16 complex with 

1053 pCp-Cy5-labelled RNA sequences of increasing length measured by 

1054 fluorescence polarization (Figure 3D). 

1055 

1056 Table 2: K d , k on and k off values of nsp10-nsp16 interaction measured by 

1057 Octet biolayer interferometry (Figure 6A) in the presence or not of SAM or 

1058 SAH. 

1059 

1060 Table 3: Measurement of nsp10/nsp16 2-O-MTase inhibition (IC S o) by 

1061 sinefungin, SAH and cap analogues, calculated from the dose response 

1062 curve of Figure 7B. 

1063 

1064 Table 4: K d values of nsp10/nsp16 and 7m GpppRNA interaction in the 

1065 presence or not of SAH or sinefungin measured by fluorescence 

1066 polarization (Figure 7D). 
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Table 1: K d values of the interaction of the nsp10/nsp16 complex with pCp-Cy5- 
labelled RNA sequences of increasing length measured by fluorescence 
polarization (Figure 3D). 


RNA 

K d (|JM) 

7m GpppG 

>20 

7m GpppGAUU 

1.07 ± 0.09 

7m GpppGAUUUAA 

0.42 ±0.01 

7m GpppGAUUUAAGUG 

0.251 ±0.02 

7m GpppGAUUUAAGUGAAUA 

0.278 ±0.05 
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Table 2: K d , K on and K ofi values of nsp10-nsp16 interaction measured by Octet 
biolayer interferometry (Figure 6A) in the presence or not of SAM or SAH. 


Condition 

K d interaction 

nsp10/nsp16 

(pM) 

Kon (M V) 

K 0 ff (s 1 ) 

Without SAM, SAH 

1.24 ±0.015 

3.77x10 3 ± 4.35 

4.68x10' 3 ± 0.016x1 O' 3 

With 100 [jM SAM 

0.15 ±0.001 

3.98x10 3 ± 2 

6.23x10" 4 ± 0.031x10" 4 

With 20 |jM SAH 

0.48 ±0.003 

3.10x10 3 ±1 .97 

1 ,51x10‘ 3 ± 0.004x10‘ 3 

With 100 mM SAH 

0.19 ±0.0013 

3.3 x10 3 ± 2.2 

6.51x10" 4 ± 0.041x10" 4 
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Table 3: Measurement of nsplO/nspI6 2-O-MTase inhibition (IC 5 o) by 
sinefungin, SAH and cap analogues, calculated from the dose response curve 
of Figure 7B. 


Inhibitor 

IC 50 (MM) 

sinefungin 

7.4 ±0.9 

SAH 

7 ± 0.4 

GpppG 

274.4 ± 140 

7m GpppG 

45 ±3.5 
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Table 
or not 


4: K d values of nsp10/nsp16 and 7m GpppRNA interaction in the presence 
of SAH or sinefungin measured by fluorescence polarization (Figure 7D). 


Condition 

K d nsp10/nsp16/7mGpppRNA 


(pM) 

Buffer alone 

0.47 ±0.03 

With 20 (jM SAH 

0.11 ± 0.01 

With 20 |jM sinefungin 

0.09 ±0.009 
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